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Abstract. Relativistic Keplerian orbital frequency (vk) and related epicyclic frequencies (radial v, , vertical vg) play an important 
role in physics of accretion discs orbiting Kerr black holes and can by resonant or trapping effects explain quasiperiodic 
oscillations observed in microquasars. Because of growing theoretical evidence on possible existence of naked singularities, we 
discuss behaviour of the fundamenal orbital frequencies for Keplerian motion in the field of Kerr naked singularities, primarily 
in order to find phenomena that could observationally distinguish a hypothetical naked singularity from black holes. Some 
astrophysically important consequences are sketched, namely the existence of strong resonant frequency for all Kerr naked 
singularities, with radial and vertical epicyclic frequencies being equal and given by the relation o» sr = cf 2 Va 2 - 1 (a 2 + l)" 1 . 
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1. Introduction 

Quasiperiodic oscillations (QPOs) of X-ray brightness had 
been observed at low-(Hz) and high-(kHz) frequencies in some 
low-mass X-ray binaries containing neutron stars or black 



holes - for a review articles see, e.g., McClintock & Remillarc 
in the case of black hole binaries and van der Klis 
I I2000I) in the case of neutron star binaries. Since the peaks 
of high frequencies are close to the orbital frequency of 
the marginally stable circular orbit representing the in- 
ner edge of Keplerian discs orbiting black holes (or neu- 
tron stars), the strong gravity ef fects must be relevant i n 
explaining high frequency QPOs jAbramowicz et all I2004I) . 
Explanation o f QPO s in the context of discs oscilla tions 
JOkazaki et all 1 19871 iNowak & Wagonetl Il99li 1 1992b has 
been c onsidered for both the warped disc s (trapped) oscil- 
lations jKato & Fukuel Il980t iKatoj l2004ah and resonant os- 
cillati ons jAbramwvicz&KluziriaM l200lHAbramowicz et all 
l2004l). 



Usually, in the kHz QPOs the power spectrum shows 
twin peaks with frequencies correlated with their X-ray inten- 
sity and the peak separation being almost constant (see, e.g 



Strohmaver et all 



1995 iFord et all 1 19971 IZhang et all Il98 



van der Klis et allll997t) . In the case of microquasars contain- 
ing stellar mass black holes, the observed ratio of the twin 
peak frequencies is exactly, or almost exactly, 3:2 - there- 
fore some resonant effects are prob ably involved for oscillat- 
ing accretion discs in microquasars jKluzniak & Abramowiczl 



l200fi iKluzniak & Abramowiczl l200ll) '. It was shown that 
the parametric resonance of the vertical and radial oscilla- 
tions at epicyclic frequencies related to the Keplerian motion 
could be the most pro bable explanation o f the observed mi- 
croquasars phenomena jTorok et alll2005l) . On the other hand, 
the forced resonance of the epicyclic frequencies, or some 
other kind of resonance with ratios given by small integral 
numbers, e.g., 2:1, 3:1, 5:2, etc. could also explain observed 
QPOs frequencies (with the same 3:2 ratio) giving them by 
combinational ("beat") frequencies jAbram owicz & KluzniakL 
l200lHToroket alll2005HAschenbachLl2004 . Really, the puz- 
zle of this 3:2 ratio kHz frequencies is still not solved defi- 
nitely and other can didates to thei r explanation, like warped- 
disc oscillations (seelKatol l2004bh or simple p-mode oscilla- 
tions (Rezzolla, 2004), are not exluded. 



Anyway, strong observational evidence supports the as- 
trophysical relevance of disc oscillation concept introduced 
during eighties (for a review of this concept see, e.g., 
iKato et all Il998h : the mechanisms of triggering the oscilla- 
tions in epicyclic frequencies were treated succefully both 
for thin (see, e.g.. IKato e t all Il9 98|) and t hick discs (e.g., 
Matsumoto et all fl989t lAbramowicz et all 120031 iRezzollal 



2004). Nevertheless, sophisticated three dimensional magne- 
tohydrodynamic simulations (3-MHD) of accretion flows, usu- 
ally does not show any twin peak kHz QPOs resembling those 



1 Interestingly, the same 3:2 ratio seems to be present also in the 
case of neutron stars sources, in dicating the same o rigin of the ob- 
serv ed quasiperiod ic oscillations jAbramowicz et all 2003; see how- 
ever, iBelionfetaiJ 2004). 
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observed Jlgumenshchev et all l2003t iDe Villiers et all 120031 
and others) - only very recentlv lKatol d2004cl) reports some 
view of the 3:2 twin peaks in 3-MHD s imulations. In addi tion, 
it has been shown, also very recently, bv lBursa et aljj2004l) that 
the possible resonant oscillations of the torus may be directly 
observable in X-ray modulation when occur in the inner parts 
of accretion flow around black hole or neutron star (even if the 
source of radiation is steady and perfectly axisymmetric). 

Apparently, the vertical and radial epicyclic frequencies of 
the Keplerian motion play a crucial role for both thin, Keplerian 
discs and thick, toroidal discs. Their properties were exten- 
sively studied in the case of accretion discs orbiting Kerr black 
holes in works mentioned above and in many others, and re- 
main to be very hot outstanding topic in astrophysics of recent 
days. On the other hand, it is natural to extend the concept of 
disc oscillations in the epicyclic frequencies around other phys- 
ical objects. 

According to the cosmic censorship hypot hesis dPenrosel 

and the uniqueness black-hole theorems kartell 1 1973ft . 
the result of the gravitational collapse of a sufficiently mas- 
sive rotating body is a rotating, Kerr black hole, rather than 
a Kerr naked singularity. Altough the cosmic censorship is a 
plausible hypothesis, there is some evidence against it. Naked 
singularities arise in various mode l s of spherically symmet- 
ric co llapse (e.g, lLake & Zanniasl 1 199(1 I.Toshi et all 120021 
12004 . In modelling the collapse of rotating stars, it had been 
shown that in some situations the mass shedding and gravi- 
tational radiation will not reduce the angular momentum of 
the star enough to lead to the formation of a Kerr black hole 
JlVTiller & de FeliceL Il984 . Candidates for the formation of 
Kerr naked singularities with a ring s ingularity were found in 
scenario of ICharlton & Clarke] Jl990l) and some 2D numerical 
models of collapsing, rotating supermassive objects imply that 
a Kerr-like naked si ngularity could develop from objects rotat- 
ing rapidly enough JNakamura et al.lTl987l) . 

It is generaly believed that black holes are stable 



singularities (Bardeenl Il973l iThornel Il974l IWald 


Il97j 


ICohen & Gautreaulll979l 


Israeli 


1986; 


de Felice & Yu 


. 1986). 



However, recently presented gedanken experiments concern- 
ing electrically charged, Reissner-Nordstrom black holes put 
some doubts on this believe. It was shown that a charged 
test particle radially falling into a nearly Reissner-Nordstrom 
black hole may tran sfer it into a Reiss ner-Nordstrom naked 
sing ularity jHubenvl 1999; see however, Ouinn & Waldl 1999 
for a more detailed analysis). Further, it was shown that an 
extreme Reissner-Nordstrom black hole could be transfered 
into a Kerr-Newman naked singularity by capturing a flat 
and electrically neutral spinning body that pl unges in radially 
with its spin aligned to the radial direction Jde Felice & Yul 
l200lh . Moreover, possible existence of naked singularities is 
supported by general mathematical studies concerning scalar 
field s around Reissner-No rdstrom naked singularities (see, 
e.g., IStalker & Shadll2004 . 

Because the cosmic censorship hypothesis is far from be- 
ing proved, naked-singularity spacetimes related to the black- 
hole spacetimes with a nonzero charge and/or rotation param- 
eter could be considered conceivable models for some exotic 



Galactic binary systems or, on much higher scale, of quasars 
and active galactic nuclei, and deserve some attention. Of par- 
ticular interest are those effects that could observationally dis- 
tinguish a naked singularity from black holes. Therefore, we 
shall discuss here in detail the properties of the vertical and 
radial epicyclic frequencies of the Keplerian circular motion 
in the field of Kerr naked singularities in order to find as- 
trophysically relevant differencies between the black-hole and 
naked-singularity cases. Indeed, we will show that in the Kerr 
naked singularity spacetimes a new, and astrophysically very 
important effect exist, which does not occur in the case of 
black holes. Namely, for any Kerr naked singularity a sta- 
ble circular geodesies exist on which the vertical and radial 
epicyclic frequencies are equal indicating thus possibility of 
an extremely strong resonance and instability of the Keplerian 
accretion discs. 

In Sect. [2] we present the well known formulae giving 
the vertical and radial epicyclic frequencies of the Keplerian 
circular orbits in the Kerr spacetimes. In Sect. [3] we discuss 
properties of these frequencies for both black-hole and naked- 
singularity spacetimes. We focus attention on the extrema of 
both the radial and vertical frequencies, if they exist. In sect.^J 
the naked-singularity case is discussed in detail and the strong 
resonant frequency, when the epicyclic frequencies are equal, 
is given as a function of the naked singularity rotational param- 
eter. In Sect.[5]we present some concluding remarks. 



2. Epicyclic oscillations of Keplerian discs 

In the case of oscillating Keplerian disc three orbital frequen- 
cies are relevant: Keplerian orbital frequency vk = Qk/2tt, 
radial epicyclic frequency v r = <w,./27r, and vertical epicyclic 
frequency vg = u)g/2jT. For discs orbiting Kerr black holes or 
naked singularities corresponding angular velo cities Qk, cor, cog 
are gi ven 2 by the well known formulae (e.g.. lNowak & Lehn 

E22I. 



£2k 



' GM S 

V r G ' 



1/2 



+ a)' 



2 r.2 

where 

a r {x,a) = 1 - 6x~ l + 8ax~ 3/2 - 3a 2 x~ 2 , 
a g (x,a) = 1 - 4ax~ 3/2 + 3a 2 x~ 2 , 

and x is introduced as dimensionless radial coordinate 

x = r/M. 



(1) 

(2) 
(3) 



(4) 
(5) 



(6) 



2 We use Boyer-Lindquist coordinates, (t, r, 9, 0). We rescale the 
central object mass with M = GMq/c 2 = ro and the central object 
angular momentum with a = 7qc/GM 2 . Here, the parameters M and 
Jo give the mass and the internal angular momentum of the Kerr black 
hole or naked singularity. 
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In the limit of the Schwarzschild black holes (a = 0), we arrive 
at 



a T (x) = 1 - -, 
x 

a g (x) = 1 , 

so that O k (x) = u>g (x) . 

In the field of Kerr black holes (a + 0), there is 



(7) 
(8) 



Qk (x, a) > a>e (x, a) > a> r (x, a) (9) 

in the range where the frequencies are well defined (Fig.^- 
right panel). 

The properties of Q. k, u> r , u>g for Ker r black-hole spacetimes 
are rewieved, e.g., in lKato et al 1 dl998l) . We can summarize that 

- the Keplerian frequency is monotonically decreasing func- 
tion of radius for all the range of black hole rotational pa- 
rameter a e (-1, 1) 3 in astrophysically relevant radii 
above the photon orbit (Fig.^- left panel) 

- for slowly rotating black holes the vertical epicyclic fre- 
quency is monotonically decreasing function of radius in 
the same radial range as well, however, for rapidly rotat- 
ing black holes this function has a maximum (Fig.[2- right 
panel) 

- the radial epicyclic frequency has a local maximum for all 
a € (-1, 1) (Fig.[l]- right panel) 

For Kerr naked singularities the behaviour of the epicyclic 
frequencies is different. We show in the next sections that the 
vertical frequency can has two local extrema and the radial one 
even three local extrema. For completeness, we shall give de- 
tailed discussion of the properties of the functions Qk(X a), 
oj x {x, a), and u>g(x, a) for both naked singularities and black 
holes. 

Obviously, all three frequencies have the general 

form, 

(GM \ m 



f (x, a) 



•r, J 



32.3| — | f{x, a) kHz. 
\M e ) J 



(10) 



3 Here and henceforth values of a > correspond to corotating or- 
bits, while a < give counterrotating orbits. 



Thus for the reader's convenience we express the frequency as 
v [Hz] M/(1OM ) in every quantitative plot of frequency de- 
pendence on radial coordinate (|6j, i.e., displayed value is the 
frequency relevant for central object with mass of 10M o , which 
could be simply rescaled for another mass just dividing the dis- 
played value by the respective mass in units of ten solar mass. 

3. Properties of the epicyclic frequencies 

First, it is important to find the range of relevance of the func- 
tions Q K (.i, a), oj r (x, a), and u>g(x, a) above the event horizon 
located at 



x + = 1 + Vl -a 2 (11) 

for black holes, and above the ring singularity located at 

x = (6 = n/2) (12) 

for naked singularities. 

The circ ular geodesies in the fi eld of Kerr black holes were 
discussed in iBardeen et alJ dl972l) . while in the case of Kerr 
naked si ngularities the circular geodesies were discussed in 
LStuchlfkl(ll980l) . 

We can summarize that circular geodesies can exist in the range 
of 



x e (xph (a), °o) , 
where 



x ph (a) = 2 



s | — arccos(-a) 



(13) 



(14) 



gives loci of photon circular geodesies. Stable circular 
geodesies, relevant for the Keplerian, thin accretion discs ex- 
ist in the range of 



x e (x ms (a), oo) , 



(15) 



where x ms (a) denotes radius of the marginally stable orbit, de- 
termined (in an implicit form) by the relation 



1 - 6x 1 +B,ax- y2 -3a 2 x- 2 = 0, 



(16) 
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which coincides with the condition 



a r (x, a) — 0. 



(17) 



For toroidal, thick accretion discs the unstable circular 
geodesies can be relevant in the range 

*mb < Xin < x < x ms , (18) 
being stabilized by pressure gradients in the tori. Here, 

(19) 



Xmb — 2 - a + 2 VT~ 



is the radius of the marginally bound circular geodesic 
that is the lower limit for the inner edge o f thick discs 
llKozlowski et~aflll978llKrolik & Hawlevll2002h . 

Clearly, the Keplerian orbital frequency is well defined up 
to x = x p h(a). However, cj r is well defined, if a r >0, i.e., at 
x>x ms (a), and w r (x) = at x ms . We can also show that for 
x>x p h, there is ag>0, i.e., the vertical frequency a>g is well 
defined at x > x p h. 

3.1. Local extrema of epicyclic frequencies 

Denoting by fl^, % r , %g the local extrema of Keplerian vk and 
epicyclic v r , vg frequencies, we can give the extrema by the 
condition 



d d 

— Vi =0o — V; = 

or ox 



for??;, ie{K,r,<9), 



(20) 



where x is dimensionless coordinate @. From we find 

that the corresponding derivatives 4 are 



K 2 



3 IGMu y/x 



re' (x^ + o) 1 



2A 



3 

where j 6 \r, 9} 



(x 3 / 2 + a) 



and 



1 a a 



a a- 



(21) 

(22) 
(23) 

(24) 
(25) 



Clearly, Q.^<0 for x>0, i.e., the Keplerian frequency is a 
monotonically decreasing function of the radial coordinate for 
any value of the rotational parameter a. 

Relations i2Q\ and ( 1221 imply the condition determining ex- 
trema fij (a) of the epicyclic frequencies: 



A (x, a) = - 



1 



2x 3 / 2 



+ a 



■ Qfj (x, a) j G {r, 8). 



(26) 



Because we have checked that in the case of counterrotating 
orbits (a < 0) the extrema 'Rg are located under the photon cir- 
cular orbit and the extrema !R r are just extension of the % x for 
corotating case (a < 0), we are focused mainly to the case of 

4 After introducing ' as d/dr. 



corotating orbits in the next discussion. In Figs 12® we show 
curves <R r k (a), k € {1,2,3} (ft/ (a), 1 6 {1,2}) implicitly de- 



termined by the relations \26l 
branches of the solution of (I26i 



indexes k, 1 denote different 



The radial epicyclic frequency has one local maximum for Kerr 
black holes 



- 1 < a < 1, 



(27) 



but it has two local maxima 5 and one local minimum for Kerr 
naked singularities with 



1 < a < a c (r) — 1.025, 

and again one local maximum for 



a > fl c( ,.) and a < -1. 



(28) 



(29) 



The vertical epicyclic frequency has a local maximum at x > x P h 
for Kerr black holes with 



a > a 



ph(fl) 



= 0.748, 



and at x > x ms for 



a > a 



ms (fl) 



= 0.952. 



(30) 



(31) 



The local maximum of oj b (x, a) is relevant in resonant effects 



for a > a 



Note that Kg has a maximum at 



^max CRu) 



0.852, 



(32) 



therefore, the situation with function a>g (x, a) is more compli- 
cated than it seems to be indicated by Fig. ^ f° r high values 
of the black hole rotational parameter a, curves cjg (x, a) cross 
each other as is shown in left panel of Fig.|4]- Fig.^does not 
show such detail because of hi-spacing between curves. 
In the Kerr naked singularity spacetimes, the function tog (x, a) 
has a local minimum and a local maximum for 



1 < a < a Cs = 1.089, 

and it has no astrophysically relevant local extrema for 



a > a Cfl and 



a < -1. 



(33) 



(34) 



Using properties of ft r (a) and %g (a), we can conclude that 
there exist two qualitatively different types of behaviour of the 
epicyclic frequencies in the Kerr black-hole spacetimes and 
three qualitatively different types of their behaviour in the Kerr 
naked-singularity spacetimes. Examples of the behaviour of the 
epicyclic frequencies for Kerr black holes are given in Fig. ^ 
(see also Fig. lA.U . 

An example of the behaviour of the epicyclic frequencies 
in Kerr naked-singularity spacetimes is shown in Fig.0](right 
panel) for the case when all the local extrema mentioned above 
are present, while for an example of the case when the num- 
ber of the local extrema is lowest see Fig. |5] The complete 
set of figures representing systematically the evolution of the 
character of the epicyclic frequencies with rotational param- 
eter growing is included in the Appendix which consists of 
Figs IA.1I (black holes) and IA.2I (naked singularities) and fold 



We distinguish which is the global one in Sect.|3 
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Fig. 2. The locations ft^ of the radial epicyclic frequency local extrema. On right panel detailed view is shown. [Here, in the next 
Fig. and hencefort we use the following convention for both kind of extrema of the radial (9?/) and the vertical CR g ) epicyclic 
frequencies: odd or missing superscript denotes a local maximum and even-numbered one means a local minimum]. For the 
answer to question whether in the case of naked singularities the maximum is global one see left panel in Fig. (HI 
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Fig. 3. The locations %i of the vertical epicyclic frequency local extrema. The right panel gives an exact information about 
position of the important points. 



also the evolution of derivatives d22t and of the ratio vg/v r of 
the epicyclic frequencies. This set of figures represents classi- 
fication of the Kerr spacetimes according to the properties of 
the epicyclic frequencies which is fully given in the Sect. [5] 
Note that in the black-hole case it is important to distinguish 
the cases when the local maximum of vq (x, a) is located above 
x ms , and under x ms . 

Clearly, the behaviour of the epicyclic frequencies substan- 
tially differs for Kerr naked singularities in comparison with 
Kerr black holes. We discuss some consequences of this differ- 
ent behaviour for naked singularities in Sect.|4] 

Now, we focus our attention on some properties of the 
epicyclic frequencies that are very important for treating the 
resonant oscillation phenomena from the observational point 
of view. 



3.2. Ratio of epicyclic frequencies 

The ratio of epicyclic frequencies vg and v, is important 



in a well defined way 


for s 


omei 


nodels 


lAbramowicz et al.[l2004llKatol 


l2004aj 


It is well known (see, e.g., 


Kat 


a et al.l 


1998) 



black holes (-1 < a < 1) the inequality 
o» r (x, a) < a>g (x, a) 
holds, i.e., the equation 

w r (x, a) = tag (x, a) 



(35) 



(36) 



does not have any real solution in all range of black hole rota- 
tional parameter a e (-1, 1) and 



Vg 



> 1 



(37) 



for any Kerr black hole. Furthermore, this ratio is a monotonic 
function of radius for any fixed a e (-1, 1) (see Fig. [5]- left 
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Fig. 4. Left panel illustrate "unlikely" effect resulting from the existence maxima of Kg (the point E on Fig. [3}- Curves vg{r) 
after a = 0.852 cross each other (curves differ in rotational parameter here by 0.05), see also Fig^for comparison. Right panel 
displays an example of epicyclic frequency behaviour for Kerr naked singularity with a = 1 .009; all allowed extrema are present. 
Note that the minimum ve is very close but not identical with the point of contact (which is also present - see Sect. l4.2l for details). 



panel). However, the situation is different for Kerr naked sin- 
gularities, see Sect. 14.21 

3.3. Implications for the orbital resonance models in 
the field of Kerr black holes 

The orbital resonance m odels for QPOs proposed by 
Abramowicz & Kluznia k JAbramowicz & Kluzniakl l200lt 
lAbramowicz et all 120041) are particularly based on resonance 
between epicyclic frequencies which are excited at a well de- 
fined resonance radius r p:q given by the condition 

^ (a, r, q ) = I, (38) 
Oh v ' q 

where p:q is 3:2 in the case of parametric resonance, and 
arbitrary rational ratio of two small integral numbers (1, 2, 
3. . . ) in the case of forced resonances. Notice that in the 
case of arbitrary forced resonance also the combinational 
("b eat") frequencies could be observed including the 3:2 ra- 
tio ( lAbramowicz & Kluzniakl l200lHTorok et all l2005l) . Such 
resonance radii are monotonically decreasing functions of the 
rotational parametr a (see Fig. [5] - left panel). Resulting res- 
onant frequencies are given gener ally as a linear com bination 
of epicyclic frequencies at r p:q . In iTorok et al.l J2005I) it is re- 
ported that the resonant frequencies (both observed frequen- 
cies, the upper and the lower) are not monotonic function of 
a for the case of 3:1 and 5:1 (5:2) forced resonance models, 
w hile for other resona nce model ratios (e.g, 2:1, 3:2) discussed 
in ITorok et alJ J2005t> . it is a monotonic function of the rota- 
tional parameter a. 

For the observational consequences, it is important to de- 
termine the limiting value of the frequency ratios ojg/cj r - pj q, 
which separates the monotonic and nonmonotonic dependence 
of the resonant frequencies on the rotational parameter. 

Indeed, this monotonicity of some resonant frequencies is 
resulting from nonmonotonic character of the epicyclic fre- 
quencies. It is known that the radial epicyclic frequency has 



a local maximum at r max ( r ) = "R T for a e (-1, 1) and its value 
Vr(max) (a) is increasing with the rotational parameter (see Fig.^ 
- right panel). Moreover, outside its maxima it is monotonically 
decreasing with the radius. From the left panel of Fig.[5]we con- 
clude that r p:q (a) must be a monotonicaly decreasing function 
of a [if the horizontal line representing some ratio p:q is fixed, 
then this figure implies necessarily monotonically decreasing 
function r p -„(aj]- Because of this, the resulting resonant fre- 
quency (which is just multiple of the radial frequency) must 
be monotonically increasing for r p:q located outside (or at) the 
maximum of the radial epicyclic frequency. For Schwarzschild 
black holes the ratio between the epicyclic frequencies at the 
radius of maximal radial frequency is exactly vg/v r = 2 (x = 8) 
and then it is slighty changing with the rotational parameter 
growing, reaching the value vgjv x ~ 1.8 for extremely rotating 
Kerr black hole (a - I). This gives the limit in the sense that for 
p:q> 1.8, the radius r p . q is surely located above the locations 
of the maximum of v r . 

On the other hand, an analogical consideration shows that vg 
(or v r , if the resonance condition fl38i is satisfied) is surely de- 
creasing with the rotational parameter, if r p:q is located under 
the location of the maximum of vg. The ratio of the epicyclic 
frequencies at the maximum of vg is shown on the middle panel 
of Fig. Its minimum is reached for extremely rotating Kerr 
black holes at vg/v r (r max ) ~ 2.18. It means that for black holes 
the resonance is surely nonmonotonic if p/q > 2.18. 

It is clear from discussion presented above that the limit 
for nonmonotonicity must be located between the values of 
p/q£ (1.8, 2.18). We checked numerically the loci of eventu- 
all nonmonotonicity and find that the limit is very close to the 
upper value, i.e., nonmonotonicity of the function vg/v r (r p:q , a) 
in dependence on a is relevant for forced resonances with 

p:q>2.U. (39) 

Figure |5] (right panel) illustrates this limit by examples of 
behaviour of vg (a, r p:q ) for three different forced resonances 
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Fig. 5. Left panel shows the behaviour of ratio vg / v r of the epicyclic frequencies - curves are spaced by 0. 2 in rotational parameter. 
Middle panel gives the information about ratio between epicyclic frequencies at maxima of vg, while right panel illustrates the 
examples of behaviour of the frequency vg for three cases of forced resonances with p/q = 5, 3, and 2.2. 



which embody the nonmonotonicy (in the sense described 
above). 

4. Some important implications for the 

discoseismology of Kerr naked singularities 

Without going into details, we expose on simple examples that 
physics of the wave propagation and oscillations in thin accre- 
tion discs around the Kerr naked singularities is considerably 
different in comparison with the Kerr black hole case. 

4.1. Modified propagation of inertial-acoustic and 
inertial-gravity waves 

The systematic classification of oscillations and waves in ac- 
cretion discs according t o the restoring fo rces acting on them 
was summarized, e.g., in lKato et alJ Jl99^) . 

They introduced the following modes for disc oscillations: 

- The higher frequency modes, i.e., inertial-acoustic waves, 
so called /9-modes. 

- The lower frequency modes, i.e., inertial-gravity waves, so 
called g-modes. 

- The inertial-acoustic modes with no node in vertical direc- 
tion are particularly called the corrugation waves. 

We remi nd that the wave propagation is given by the dispersion 
relation (bkazaki et al.lll987tlKatolll989tlNowak & Wagoned. 
ll99lUl992h " 

a) 2 - n Yd 2 , - = a? k 2 c 2 , (40) 

where co is the angular velocity corresponding to the frequency 
of propagated wave, T is a correction factor of the order of 
unity, k r is the wave vector and c s is the vertically averaged 
sound speed. For the case of axially symmetric (m = 0) funda- 
mental mode (« = 0) is this relation reduced to 

to 2 = a) 2 + k 2 c 2 . (41) 

In this simple case the only waves of p-mode (co 2 > a> 2 ) can 
propagate in the radial directio n and could be in t he black hole 
case separated into two groups dKatoet al.Lfl998i) : 



- The waves with 6 (J 2 > <j^ lax can propagate through the en- 
tire region of the disc. 

- For the waves with to 2 < <w^ ax the propagation region is sep- 
arated into two discrete portions (see Fig. [6]- left panel): 

- Inner region r < r\ . 

- Outer region r > r%, 

where r\, ri are the solutions of aj 2 -a> 2 ; n- waves with 
it) 2 < w^ ax cannot propagate between r\ and n, however they 
can propagate outside this region. Because of physical proper- 
ties of the inner edge of the disc, the n waves could be trapped 
in the inner region, i.e., they are reflected be tween the inner 
edge and r\ as shown bv lKato & Fukuel Jl98(j) . 

The overtones (n ^0) in equation (I40> allow the inertial 
acoustic p-waves as well as inertial gravity g-waves. Generally, 
p-waves in this case must have frequency or > n TQ^, while 
for g-waves the inequality to 2 < to 2 must hold. 
The g-waves which could be trapped (see Fig.|S]- right panel) 
in the region around r max(1 ) = 9? r are insensitive to changes of 
the disc structure, if the condition of moderately geometrically 
thin disc is satisfied dOkazaki et allll987llNowak & Wagoner, 
ll992HPerez eTailll997l) . 

In the case of Kerr naked singularities, the well known 
picture of the wave propagation, shortly reminded above, is 
changed due to the different behaviour of the epicyclic fre- 
quencies. As shown in Sect, fj] the radial epicyclic frequency 
has two global maxima for Kerr naked singularities with 
1 <a <a cr = 1.025. Thus the propagation of the fundamental 
axisymmetric (n = 0, m = 0) 7r-waves is changed (of course, the 
condition or >u> 2 must be satisfied again): 

- The waves with 7 a> 2 > cu 2 MAX can propagate through the 
entire region of the disc as well as in the black hole case. 

- The waves with oj 2 < w 2 MAX could be divided into three 
groups 



6 Here and hencefort in relation to the wave propagation for black 
holes a> max denotes the maximum of a> T at R r (a) (Fig.|2J in the black- 
hole range of the rotational parameter. 

7 In the relation to the wave propagation in discs around Kerr naked 
singularities, (jmax denotes the global maximum of w r , while a> max 
denotes a local (lower) one, see Fig.|H| 
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Fig. 6. rAfter lKato et~ailll998l) The simplest examples of wave propagation in black-hole thin accretion disc: for fundamental 
axisymmetric mode (n = m= 0) see left panel, while the first axisymmetric overtone (n = 1, m= 0) is shown on right one. Both figs 
are plotted for the Schwarzschild black hole. 
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Fig. 7. The same examples of wave propagation as in Fig.|6] but for naked-singularity thin accretion disc (for the simplicity in 
degenerated case of a = 1.012). Fundamental axisymmetric modes appear on left panel and the first axisymmetric overtone on 
the right one. 



- For waves with a) > w max the propagation region splits 
into two (inner and outer) regions as in the case of black 
holes. 

- For waves with w m j n <a>< a> max the propagation region 
splits into three different portions from which two could 
contain trapped waves. 

- For waves with oj m \ n > ai the propagation region splits 
again into inner and outer region as in the black hole 
case. 

This modified wave propagation is illustrated in Fig. 
(left panel) for the simple degenerated case a =1.012 when 

Wmax = WMAX- 

In Fig.Q(right panel) we illustrate that analogically to the 
situation described above for trapped 7r-waves of fundamen- 
tal mode, the propagation of trapped g-modes (« + 0) splits 
into discrete portions in the case of naked singularities with 
a < a cr = 1.025. 



For the Kerr naked singularities with rotational parameter 
a > a cr , the radial epicyclic frequency has the same properties 
as in the black hole case and the wave propagation is formally 
the same as for black holes. 



4.2. Strong resonant frequency 

It is shown in Sect. [3] that for Kerr naked singularities with 
a > a C (g) = 1 .089 the behaviour of epicyclic frequencies is for- 
mally similar as for Kerr black holes. However, for any naked 
singularity with a > 1, the epicyclic frequencies (|2]|3j can sat- 
isfy the equality condition 



a>i(a, x) = a>e(a, x) (42) 
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Fig. 8. Left panel shows the radial epicyclic frequency at particular extrema as function of rotational parameter a between a — 1 
and a = 1.025. We can see that for a < 1.012 the global maximum is situated at (a) while for a > 1.012 it is at ft? (a). Middle 
panel exposes the behaviour of function vg/v r (x) typical for Kerr naked singularities. With increasing parameter a the local 
maximum of this function is shifted to higher radii. Right panel illustrates that the value of maximum vg/v r is rapidly decreasing 
with the rotational parameter growing. Such behaviour of function vg/v r (x, a) implies the important consequence for resonance 
models: in the case of slowly rotating Kerr naked singularities an eventual resonance orbit r p: „ (with p, q being small integral 
numbers) is defined unambiguously - the upper limit for this unambiguity in the case of p:q = 3:1, 2:1 and 3:2 is denoted. 



giving strong resonant phenomenon 8 which occurs at the crit- 
ical radius 



x sr — a" (a > 1). 



(43) 



This means that for any Kerr naked singularity the epicyclic 
frequency ratio ve/v r (r) is a nonmonotonic function which 
reaches the value 1 at the point given by J43l > (see Fig. [8]- mid- 
dle panel). The loci of this point is compared with locations 
of some other important points as show in Fig. [9] - the right 
panel, while the left panel shows an example illustrating radial 
extension of the strong resonant phenomenon. 

Using the relation J43l > in (|2j|3} we find the strong resonant 
frequency which in terms of the corresponding angular velocity 
reads 



0J st = to r (a, a~) = u>o (a, a ) 



(GM 



\ G ' 



1/2 



a 2 (a 2 + lY 



and the frequency can be expressed in the form 



Vsr = 32.3 



I M ! a 2 (I 



2 (1 +a 2 ) 



kHz. 



(44) 



(45) 



We note that the strong-resonance phenomenon represents a 
crucial difference of Kerr naked singularities from the case of 
Kerr black hole for which the ratio u>g/u) r (r) is determined as 
a monotonic function for fixed a (Fig. [3]- left panel). Notice 
that for Kerr naked singularities with any rotational parameter, 
the ratio (x>gju x < (we/w r ) max being fixed can appear at three 
different radii 9 . Note that this kind of behaviour is implied by 
the existence of the strong resonance frequency (a>g = 0)r) for 
any Kerr naked singularity, i.e., it is not restricted to the cases 
when local extrema of cog, oj r exist. 

The behaviour of the epicyclic frequency ratio a>g / a> T (r) 
typical for Kerr naked singularities is shown in Fig. [8] (mid- 
dle panel). In right panel in Fig.|S]we plot the value of the local 



l43l 



Intuitively clear attribution is well founded in the last subsection 



9 This can be important in discussion of observational aspects of 
resonant phenomena. 



extrema of the ratio u)g/u> T as a function of the rotational param- 
eter a. We see that for high values of the rotational parameter, 
the radial and vertical epicyclic frequencies are very close each 
other in large radial range around r sr . The example is given in 
Fig. El(left panel). 

We plot the strong-resonance frequency as a function of the ro- 
tational parameter in Fig.^|(left panel.) The strong-resonance 
frequency is approaching zero value for extremely rotating 
Kerr black hole and has a maximum for naked singularities 
having rotational parameter 

a src = 1.207, (46) 

with the corresponding value of the epicyclic frequency being 
determined by the relation 

v src = 6.1 kHz. (47) 

In the middle panel of Fig.[H)| location of the critical radius 
x sr = a 2 is shown together with location of the marginally stable 
orbit and extension of the instability 10 region of r where the 
difference between values of the radial and vertical epicyclic 
frequencies is smaller then 1%. However, for the values of ro- 
tational parameter a — 1 , it is more convenient to express the 
region of the disc with 1 % difference of the epicyclic frequen- 
cies in terms of the proper radial distance f, which has direct 
physical meaning. There is 



r 

•J l „ 



(48) 



where g rr denotes the radial metric coefficient of the Kerr met- 
ric in the standard Boyer-Lindquist coordinates; the distance 
is measured from the inner edge of the thin discs, located at 
r ms . The result is represented by right panel in Fig.^3 we nave 
found that location of the strong-resonance frequency is clos- 
est to the inner edge of the Keplerian disc for naked singularity 
with 



fl f = 1.105. 



(49) 

10 The meaning of the term instability used here is discussed in the 
last subsection l4.3l 
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Fig. 9. Left panel shows the behaviour of epicyclic frequencies for a = 2.3, the region where frequencies are identical with 
accuracy of 1 % is denoted. Right panel illustrates the location of the strong resonant frequency (dotted curve a 2 ) in relation to 
extrema of epicyclic frequencies (thin curve thick one The critical radius is always located (in radial order) between the 
first maxima of a> r and the minima of tog (if these exist). The notation of the important points [A. . . J] is in accord with Figs.|2l[3] 



For this value of the rotational parameter the strong resonant 
frequency is 

v sri n=4.3^J kHz, (50) 

what is about 70% of maximum at fl srmas given by J47I . 
However, we remark that the critical radius is always located 
outside from the innermost part of the disc (see Fig. [9] right 
panel). 

4.3. Possible instability of the accretion disc around 
Kerr naked singularities 

The orbital resonance model feluzniak & Abramowiczl 120021) 
demonstrates that fluid accretion flows admit two linear quasi- 
incompressible modes of oscillations, vertical and radial, with 
corresponding eigenfrequencies equal to vertical and radial 
epicyclic frequencies for free particles. In a particular model 
of slender torus, the general properties of this modes can be 
shown: the vertical mode corresponds to a periodic displace- 
ment in which the whole torus moves as a rigid body up and 
down the equatorial plane and each fluid element has a vertical 
velocity that periodically changes in time, but does not depend 
on the position. The frequency of the vertical mode is equal 
to the vertical epicyclic frequency that a ficticious free particle 
orbiting at the circle of maximum pressure in the torus equilib- 
rium position would have. Behaviour of the radial mode is sim- 
ilar to the vertical one, a nd in the linear regime these two modes 
are formally uncoupled. lKluzniak & Abramowiczl J2002I) argue 
that in the case of more realistic description which includes 
non-linear effects given by pressure and dissipation, non-linear 
effects couple the two epicyclic modes what may results in a 
resonance. 

One of such possible resonances, the parametric reso- 
nance, seems to be most probable explanation of the 3:2 dou- 
ble peak kHz QPOs observed in some galactic microquasars 



jAbramowicz et alll2004l). The effect itself i s described by the 
Mathieu equation JLandau & Lifshitzll 19761) 11 

66 + u 2 g [ 1 + h cos(w,.f)] 66 = 0, (51) 

that can be formally derived by considering small deviations of 
fluid streamli nes from planar circular motio n governed by a set 
of equations jRebuscol l200llHorakLE004l) 

6r + co 2 r 6r = u 2 f r {6r, 66, 6r, 66), 

66 + u]6d = co 2 g f g (6r, 66, 6r, 66), (52) 

for a particular choice of f r and fg, corresponding to 

6f + aj 2 .6r = 0, 5d + (jj 2 e 56 = -(jj 2 e 666r. (53) 

From the theory of the Mathieu equation it is known that the 
parametric resonance is then excited when 12 

a> r v r 2 

— = — = -, n = 1,2,3.... (54) 

OJg Vg n 

The effect is strongest fo r the smallest possible value of n 
(Landau & Lifshit7t Tl97fil). Because in the field of black holes 
v,. < vg (see Sect. [3}, the smallest possible value for resonance 
is n = 3, i.e, 2 vg — 3 v,-, wh at explains very well the 3:2 ra- 
tio observed in mi croquasars dXb ramowicz &Kluzniakll2004l 
lToroketal.Ll2005l) . 

As we show in previous text, the point with radial and ver- 
tical epicyclic frequency being equal exists for any Kerr naked 
singularity. Obviously, at such point the equality (I54> is sat- 
isfied (v r /vg = 1/1 =2/«; n = 2), and the parametric resonance 
(between the radial and vertical epicyclic frequency) eventually 
excited at this point is the strongest possible parametric res- 
onance excited between the epicyclic frequencies in the field 
of Kerr naked singularities. Such 2:2 resonance must be also 

11 After denoting the time derivativ d/df by dot. 

12 We note that the same condition holds fo r internal reso nance, 
which describes systems with conserved energy <Horakll2004l) . 
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Fig. 10. The behaviour of strong resonant frequency is shown on left panel. From the middle panel we can see the evolution of 
area (around critical radius) where epicyclic frequencies are close, while right panel show the same in the proper distance to 
marginally stable orbit. 



stron ger then the 3:2 parametr ic resonance in the black hole 
case ( Landau & LifshitzLll973) . 

From this, and from the fact that the radial region with the 
epicyclic frequencies being nearly equal is rather large, one can 
expect that at this region both radial and vertical oscillations 
could be strongly amplified, leading to an instability of the ac- 
cretion disc 13 . 



5. Conclusions 

For counterrotating Keplerian orbits, properties of the epicyclic 
frequencies are the same for all Kerr black holes and naked sin- 
gularities. Radial epicyclic frequency has always a local max- 
imum, while the vertical epicyclic frequency has no local ex- 
trema at x > x pn . 

On the other hand, for corotating Keplerian orbits proper- 
ties of the epicyclic frequencies strongly depend on the rota- 
tional parameter of the Kerr spacetimes. Altough, from this 
point of view, the most important difference between space- 
times with a < 1 and a > 1 is the change of inequality 

a>o(x) > w r (x) (a < 1) — * a)g(x) > cj^x) (a > 1), (55) 

we have to distinguish different possibilities according to the 
existence and relative locations of the local extrema of the 
epicyclic frequencies. 

In the case of Kerr black holes, the classification according 
to the properties of the epicyclic frequencies is given in the 
following way: 

- BH1 (Fig.lATV) < a < 0.748 

u> T (x, a) has one local maximum, u>g(x, a) has no local ex- 
trema above the photon circular orbit x pn . 

- BH2 (Fig.EU>) 0.748 < a < 0.952 

a> T (x, a) with one local maximum, u>g(x, a) with one local 
maximum at x < x ms . 

13 Such a claim is motivated by experience from known situations re- 
lated to the parametric or forced resonance in complex non-linear sys- 
tems observed in Earth physics lLandau & Lifshitzlll97d) . Examples 
of mathematically possible resonances causing damages of bridges, 
wings, etc. with no specific physica l coupling mechanism known are 
discussed in lNavfeh & Mookl 1 19791) . 



- BH3 (Fig. IXlt) 0.952 < a < 1 

a> T (x, a) with one local maximum, U)g(x, a) with one local 
maximum at x > x ms . 

In all the cases, the function (u>g/u> T ) (x, a) which is relevant 
for resonant effects, has a monotonic (descending) character 
(Fig. IA.lt . Therefore, for a given rotational parameter a, there 
is only one radius allowed for any p : q resonance. However, 
we have shown that the resulting resonant frequencies are non- 
monotonic functions of a for p : q > 2.18, what could embar- 
rass eventual spin estimate given by some resonance models. 
In addition, the curve 'Rgia) has a local maximum at a = 0.852 
and we can conlude that the functions u>g(x, «i), a>g{x, ai) with 
a\, ci2 being fixed and higher then a 0.85 cross each other - 
this can be also of observational interest. 

In the case of Kerr naked singularities, the classification is 
given in the following way: 

- NaSl (Fig.EUO 1 < a < 1.012 

oj t (x, a) has two local maxima and one local minimum be- 
tween the maxima; 

Wr(max) (^r(in), fl) < ^i(max) (*r(out), «) < Wfl( max ), where 
X T (in) = ' x r (out) = ■ 

u>g(x, a) has one local minimum and one local maximum. 
There is x 9 ( max) = fl l e < x r(out ). 

- NaS2 (Fig.rOb) 1.012 < a < 1.024 

oj r (x, a) has two local maxima and a local minimum in be- 
tween, 

W r ( max ) (jCr(out), «) < W r ( max ) (jC r (in)) a ) < W9(m H ). 

a>g(x, a) has one local minimum and one local maximum, 
with 

- NaS3 (Fig.|A2t) 1-024 < a < 1.025 

The same as in the class NaS3, but with xg( max ) > x r ( 0u t)- 

- NaS4 1.025 <a < 1.047 

oj t (x, a) with one local maximum. a>g(x, a) with one local 
minimum and one local maximum; u>g( m . dX ) > u>r(max)- 

- NaS5 (Fig.rOli) 1.047 < a < 1.089 

The same as class NS4, but with cog( ma x) < w r ( max ). 

- NaS6 (Fig.ESfc) a > 1.089 

0) r (x, a) with one local maximum, u>g{x, a) with no local 
extrema. This class is formally similar to the class BH1, 
but with the crucial exception of the point x SI . 
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We conclude that the properties of the radial and verti- 
cal epicyclic frequencies of the Keplerian motion in the case 
of Kerr naked singularities differ substantially from the case 
of Kerr black holes and can have strong observational conse- 
quences for both resonant phenomena and stability of accretion 
discs around Kerr naked singularities. 

While the wave propagation in the case of oscillations 
in the discs around Kerr naked singularities is substantially 
different from the black hole case only for spacetimes with 
a € (1, 1.025), the strong resonance effect can occur for any 
Kerr naked singularity - the strong resonant frequency of disc 
oscillations around Kerr naked singularities arise always at the 
descending part of the function <w,(x, a) (in vicinity of the lo- 
cal minimum of a>g, if this exists), i.e., it is always located 
above the innermost part of the disc. We stress that this phe- 
nomenon represents the strongest parametric resonance be- 
tween the epicyclic frequencies possible in the field of Kerr 
naked singularities, stronger than in the case of Kerr black 
holes. Moreover, the area where the effect occurs is large what 
could finally have a strong influence on the stability of the 
Keplerian disc itself. 

It follows from the existence of the strong resonant 
frequency that the function (oJa/w r ) (x, a) has the same 
character for all Kerr naked singularities having one lo- 
cal maximum (we/w r ) max > 1 and one local minimum at 
(wfl/a^Otnin {x sl - = a 2 , a) - I which corresponds to the strong 
resonant frequency. This implies an important consequence for 
the other resonant phenomena - namely, in the range of fre- 
quencies (wfl/w r ) max (a) > GJg/cOr > 1, the resonant effects with 
the same rational ratio can occur at three different radii r p:q . 
Using the behaviour of the function {tool ^)m-dx («) (Fig. H] - 
right panel) we can conclude that three radii r p - g could occur in 
the field of Kerr naked singularities with a < 1.0012 for p : q 

< 3 : 1, a < 1.012 for p : q < 2 : 1, and a < 1.062 for p:q 

< 3 : 2. Clearly, in such situation the same resonant frequency 
ratios could be induced by very different physical phenomena 
at different parts of the accretion disc. 
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Appendix A: Classification of the Kerr spacetimes (due to epicyclic frequencies) 



a) a = 0.5 e (0,0.748) = BH1 




T U 



d) degenerated case of a = 1 




Fig. A.l. Classification of the Kerr black-hole spacetimes: behaviour of the epicyclic frequencies (left panel), their first deriva- 
tives (middle panel), and their ratio vg / v r (right panel) are shown for four representative values of rotational parameter a including 
the extreme Kerr black hole. Left margin of plots is always situated at the photon circular orbit r p h, while the marginaly stable 
orbit r ms is denoted by dashed vertical line. 
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Fig. A.2. Classification of the Kerr naked-singularity spacetimes: behaviour of the epicyclic frequencies (left panel), their first 
derivatives (middle panel), and their ratio vg/v r (right panel) are shown for five representative values of rotational parameter a, 
example of the class NaS4 which differs from the class NaS3 by the absence of the "radial pair" maximum-minimum is not 
shown. 



